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Abstract 

In the present paper we investigate incoherent p° electroproduction off com- 
plex nuclei. We derive a novel, simple expression for the incoherent elec- 
troproduction cross section in which one can clearly separate the final state 
interactions of the reaction products from the 'initial state interactions' of 
the photon that give rise to nuclear shadowing. In the special case of purely 
absorptive final state interactions we deduce from our expression the known 
Glauber result. A more realistic treatment of the final state interactions 
within a transport model is then used to compare our predictions with exper- 
imental data from the HERMES experiment. 
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I. INTRODUCTION 



High energy meson electroproduction off complex nuclei offers a promising tool to study 
the physics of hadron formation. The nuclear environment makes it possible to investi- 
gate the properties of hadrons immediately after their creation. In addition one can vary 
the energy and virtuality of the exchanged photon to examine the phenomenon of color 
transparency [1]. 

In the present paper we study exclusive incoherent p° electroproduction off 14 N and 
84 Kr in the energy regime of the HERMES experiment [2,3]. In this kinematic regime the 
coherence length of the photon's p° component is large and quantum mechanical interference 
effects play an important role. In [4] we presented a method to account for this so-called 
shadowing effect in photoproduction within a semi-classical BUU transport model which 
we now apply for electron induced reactions. We distinguish between the 'initial state 
interactions' of the photon that give rise to shadowing and the final state interactions (FSI) 
of the produced particles. The coupled channel treatment of the FSI within the transport 
model allows for a broader spectrum of FSI than the usually used Glauber theory [5]. In 
addition we account for Fermi motion and Pauli blocking during the reaction and take the 
finite lifetime of unstable reaction products into account. 

The time that the products of an elementary photon-nucleon interaction need to evolve 
to physical hadrons is called formation time. Any change in the reaction probability of a 
hadron h during its formation time affects the so-called nuclear transparency 

T A = ^ A - hX (1.1) 

if the formation length, i.e. the length that h travels during its formation time, exceeds the 
average internucleon distance (~ 1.7 fm) inside the nucleus. 

In the case of diffractive vector meson electroproduction (7*iV — > VN) the virtual photon 
initially produces a colorless gg-pair whose size is expected to decrease with increasing photon 
virtuality Q 2 [6]. As long as the gg-pair is very small, i.e. in the early stage of its evolution 
into the physical vector meson, it mainly reacts via its color dipole moment. This leads to 
a cross section that is quadratic in the gg-pair's size. At large enough energies the gg-pair 
is frozen in this small sized configuration over a distance that, because of time dilatation, 
can exceed the diameter of a nucleus. This effect leads to a large nuclear transparency for 
vector meson production at large photon energy v and virtuality Q 2 (color transparency). 
For a detailed review of color transparency see e.g. [7] and references therein. 

It is not clear whether the energy and the virtuality of the photon in the HERMES 
experiment is large enough to see an onset of color transparency. In the present paper 
we do not find evidence for a finite formation time of the produced vector meson. The 
14 N data for the transparency ratio as function of the coherence length is compatible with 
the assumption that the diffractively produced p° starts interacting with a hadronic cross 
section right after its production. Kopeliovich et al. [6] point out that this might be an 
accidental consequence of the specific correlation between the Q 2 of the photon and the 
coherence length of the p° in the HERMES data. In their work a light-cone QCD formalism 
was used to incorporate formation and coherence length effects in coherent and incoherent 
vector meson electroproduction. In addition the effect of gluon shadowing was studied, but 
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found to be negligible in the case of incoherent p production. In Ref. [6] the decrease of the 
transparency ratio due to the finiteness of the p° lifetime was discussed as well. This effect 
is automatically included within our transport model. Our coupled channel calculation 
shows in addition that the effects of side-feeding in the FSI are unimportant because of 
the kinematic cuts of the HERMES experiment. Therefore Glauber theory seems to be 
appropriate for these reactions. However, we find that for heavy nuclei the kinematic cuts 
lead to an additional reduction of the transparency ratio due to elastic FSI which scatter 
the produced particle out of the acceptance window. One must therefore be careful which 
cross section one uses to describe the FSI in Glauber theory. 

Our paper is structured in the following way: In Section II A we show how we describe 
the electron-nucleon interaction and how we account for shadowing within our model. In 
the limit of purely absorptive FSI our result for the incoherent vector meson production 
cross section turns out to be equivalent to the known Glauber result of Huefner et al. [8]. 
However, in contrast to the result of Ref. [8] our expression allows for an intuitive physical 
interpretation. The transport model itself is sketched in Section II B. The results of our 
calculations for the exclusive incoherent p° production cross section within the coupled 
channel treatment are presented in Section III in comparison with experimental data from 
the HERMES experiment. We close with a short summary and outlook in Section IV. 

II. MODEL 
A. Shadowing 

In the one photon exchange picture the interaction of an electron and a nucleon (or 
nucleus) can be reduced to the reaction of a virtual photon on the hadronic target. The 
exchange of more than one photon is suppressed by further factors of the fine structure 
constant a em . We follow the method of Friberg and Sjostrand [9] and use the event generator 
PYTHIA v6.2 [10] to describe the interaction of the (virtual) photon and a nucleon. The 
basic idea is that in a photon hadron collision the photon not necessarily interacts as a 
point particle (direct interaction) but might fluctuate into a vector meson 
(vector meson dominance) or perturbatively branch into a qq pair before the interaction 
(generalized vector meson dominance, GVMD). Within the PYTHIA model it turns out 
that the latter is very unlikely in the kinematic regime of the HERMES experiment (photon 
energy v 10 — 20 GeV, Q 2 ~ 0.5 — 5 GeV 2 ) as can be seen from Fig. 1 where we show the 
contribution of the different photon components to the total 7*V cross section at an invariant 
mass W =5 GeV. In most cases the photon will therefore scatter deep inelastically from a 
parton of the target nucleon or fluctuate into a vector meson before it reaches the nucleon. 
In the latter case the vector meson might either scatter diffractively from the nucleon or a 
hard scattering between the constituents of the vector meson and the nucleon might take 
place. The hard scattering leads to the excitation of one or two hadronic strings which finally 
fragment into hadrons. We assume that the hadrons that emerge from the edges of these 
strings, i.e. those that contain quarks not originating from the string fragmentation, can 
interact with a hadronic cross section right after the photon nucleon reaction. Consequently, 
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for the diffractively produced vector meson we do not use a formation time since none of its 
constituents arise from a string fragmentation. 

If the struck nucleon is embedded in a nucleus one has to account for its Fermi motion 
and binding energy as well as Pauli blocking of final state nucleons. In addition one has 
to be aware that the nuclear environment influences the VMD part of the photon nucleon 
interaction, since the vector meson components get modified on their way through the nuclear 
medium to the interaction point. 

We express the physical photon state (7) in terms of vector meson states \ V) and a state 
|7o) which consists of the point-like photon and the GVMD part. This makes sense if one 
assumes that the GVMD part does not get shadowed or as in our case is unimportant for 
kinematic reasons: 

i7>=(i- E ^?W>+ E f/y\n (2.1) 



The formfactor 



Fy 



Q 2 + W 2 



1 + 0.5- 



Am 2 v Q 2 



m 2 v + Q 2 ) 2 \ Wv + Q 



(2.2) 



is taken from [9] and also accounts for contributions of longitudinal photons. In Eq. (2.2) 
v denotes the energy of the photon, Q 2 its virtuality and my the mass of the vector meson 
V. W is the invariant mass of the photon nucleon system. 

For a photon with momentum ke z the wave function of the component V can be written 

as 

^\b,z)= X (by kz . (2.3) 

On the way through the nucleus the vector meson components get modified. Consider an 
ensemble of A nucleons at positions {r*j} = {(sj,Zj)} which are labeled in such a way that 
z\ < Z2 < ... < za- According to Glauber theory [11] the wave function ipy behind the first 
nucleon looks like 

^\b,z) =i> { °\b,z) -r v (b-s 1 )^\b,z 1 )e ik ^ z -^ (2.4) 
= X$) (e lkz - T v (b - sOe^V^ 2 ) 

if one neglects off-diagonal scattering from one vector meson component into another (di- 
agonal approximation). The phase factor in the second term arises from putting the vector 
meson component on its mass shell, i.e. ky = \/ v 2 — rriy and qy = k — k v is the corre- 
sponding momentum transfer. The profile function T v is related to the elastic vector meson 
nucleon scattering amplitude jy in the following way: 

Ty(b) = -Xr- [ ^ rb fy{k t )d 2 k t . (2.5) 



2%iky j 

After passing the second nucleon the wave function takes on the form 
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^\b,z) = X (b) (e lkz - F v (b- s 2 )e iqvZ2 e lkvZ - T v (b - ^(l - Y v {b - s 2 )e iqvZl e ikvZ ^ . 

(2.6) 

At position z the whole ensemble of A nucleons has led to the modification 

A A 



^ v A \b,z) = X (b) e ikz -J2 T v(b- sj)e^ J] [l - T v (b - s k )Q(z k - Zj j\ &(z - z k )e l 



ikyz 



k+3 



A A 1 

1 -J2r v (b- s,Y^ z ^ J] [l - T v (b- s k )Q(z k - Zj j\ S(z - z k ) \ $\b,z) 

3=1 k^j J 

= {l-r^(M;{r1})}^ 0) (M). (2.7) 

— * 

Hence, the physical photon state changes on its way through the nucleus to position r = (ft, z) 
according to 



17(0) = 1 - E 



2o 2 

v= p ,oj,4>,J/^ yy 



n) l7o> + £ f Fy (l - rt } (f)) |V) (2.8) 



where the nuclear profile function from (2.7) must be averaged over the positions of the 
nucleons in the nucleus 



r^(ft,^) = (o|r^(ft,^{ri})|o) 



(A), 



dzin(b, Zi) VN (1 — iay)e 



exp 



--a VN (l - ia v ) I dz k n(b,z k ) 



(2.9) 



To get to the second line of Eq. (2.9) we have followed the way of [11], i.e. making use of the 
independent particle model and the large A limit and used the optical theorem to express 
the elastic forward scattering amplitude in terms of the total VN cross section ay: 



MO) = % ^-a v {\-ia v ). 



(2.10) 



The nucleon number density is denoted by n(r) and normalized to A. 

Neglecting any influence of the FSI, the reaction amplitude for the process 7iV — > / on a 
nucleon at position f inside a nucleus changes compared to the vacuum due to these 'initial 
state interactions' of the photon: 



(/m 7 > - (f\Tm). 



(2.11) 



By comparing Eq. (2.1) and (2.8) one sees that to account for shadowing in incoherent 
photoproduction one only has to multiply the amplitude of each vector meson states \V) by 

the corresponding factor ^1 — T^\f)j. This is equivalent to the method we used in [4]. 
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We now compare our method of describing shadowing of incoherent electroproduction 
with conventional Glauber results. We therefore calculate the incoherent vector meson elec- 
troproduction cross section off a nucleus and assume purely absorptive FSI for the moment. 
After averaging over all positions r — (b, z) we get: 



<7 1 A^VA* = C-yN->VN d b dzn(b, z) 



X 



2 

1 - J dzin(b, zi)-^-(l - ioL V ) 



e iqv(zi-z) exp 



2 

~ov(l-io!y) J dz k n(b,z k ) 



x exp 



oo 

-ap el J dz'n(b, z') 



(2.12) 



In (2.12) we again neglect off-diagonal scattering, i.e. V production is triggered by the V 
component of the photon only. The last factor in (2.12) includes the inelastic VN cross 
section ay el and accounts for the FSI. 

The production process described by Eq.(2.12) can be understood as follows: The inco- 
herent vector meson production takes place on a nucleon at position {b, z). The production is 
triggered either directly by the photon or via an intermediate vector meson which was earlier 
produced on a nucleon at position (b, Zj) without excitation of the nucleus. The interactions 
of this vector meson on its way from Zi to position z are of optical potential type and leave 
the nucleus in its ground state. This propagation is described by the second exponential 
in (2.12). The interference between the direct and the indirect process leads to shadowing. 
The nucleus gets excited in the incoherent reaction at position z. The possibility that the 
vector meson is lost on its way from position z out of the nucleus is taken care of by the last 
exponential. 

Expression (2.12) formally differs only slightly from the approximate result for the inco- 
herent vector meson production cross section which is given without derivation in Ref. [11]. 
The latter can be obtained by making the following replacements in (2.12): 



iq v {zi-z) 



a 



inel 
V 



a v . 



(2.13) 



The last replacement simply means that one neglects the possibility of elastic VN scattering 
whereas the occurrence of the different phase factor is unclear. The Authors of Ref. [8] have 
given a seemingly different result for incoherent V electroproduction: 



a lA _ VA * = a~ jN ^ VN Id b< dzn(b, z)e Uv 



1 ay 



+ 2^ 



inel 



a 



V ) 



+oo +oo 

J dzin(b,zi) J dz 2 n(b,Z2) 



— OO Zl 
1 (rrinel _„cl\ 



x cos [qy( Zl - z 2 )} e -^v-<)T, 2 (b)-^ v T zl{b ) 



6 



i K) 2 

4 ^ 



+00 



} 



(2.14) 



with T 2 (6) = n(b, z')dz' and T(6) = T 2= _oo(6). Expression (2.14) can be interpreted 
as follows: The last term represents the coherent V photoproduction cross section which 
is subtracted from the inclusive V photoproduction cross section to yield the incoherent 
part. In Appendix A we show that the physically more transparent expression (2.12) is 
mathematically identical to (2.14). 



B. Transport model 

Up to now we have ignored the effects of Fermi motion, binding energies and Pauli block- 
ing of final state nucleons in the process described by (2.11). In addition we neglected the 
finite life time of the p° and took only absorptive FSI into account. These shortcomings are 
avoided when using the model of Ref. [4,12]. The propagation of the final state |/) through 
the nucleus is treated within a semi-classical transport model based on the Boltzmann- 
Uehling-Uhlenbeck (BUU) equation. The BUU equation describes the time evolution of the 
phase space density fi(r,p,t) of particles of type % that can interact via binary reactions. 
Besides the nucleons these particles involve baryonic resonances and mesons (n, r], p, K, ...) 
that are produced either in the primary reaction or during the FSI. For a particle species % 
the BUU equation can be written as: 

(0 OH OH \ 
ot + ~a? a?- WW mm = ' fM] - (2 ' 15) 

For baryons the Hamilton function H includes a mean field potential which depends on the 
particle position and momentum. The collision integral on the right hand side accounts for 
the creation and annihilation of particles of type % in a collision as well as elastic scattering 
from one position in phase space into another. For fermions Pauli blocking is taken into 
account in I co u via blocking factors. The BUU equations of each particle species i are coupled 
via the mean field and the collision integral. The resulting system of coupled differential- 
integral equations is solved via a test particle ansatz for the phase space density. For details 
of the transport model see Ref. [12]. 

The classes of FSI that are included in the transport model goes far beyond what can 
be achieved within Glauber theory. As a result the finally observed p° does not need to be 
created in the primary reaction but might be produced during the FSI via side feeding. It 
is therefore clear that a purely absorptive treatment of the FSI as in Glauber theory can 
only be used if one is sure that one has eliminated the possibility of side feeding by applying 
enough constraints on the observable (see [4] for details). 

We also stress that within our model instable particles might decay during their propa- 
gation through the nucleus. In case of the p° this means that both pions have to escape the 
nucleus without further rescattering to make the identification of the p° still possible. 
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III. RESULTS 



Before we turn to nuclear targets we first verify that the input of our model is reasonable. 
We therefore look at exclusive p° production off hydrogen. We use the same kinematical 
cuts as in the HERMES experiment [13], i.e. the final state has to consist of two oppositely 
charged pions with invariant mass between 0.6 GeV and 1 GeV. The four momentum transfer 
\t — t max \ between the virtual photon and the 7r + 7r~-pair has to be smaller than 0.4 GeV 2 
and we apply the exclusivity measure 

AE = p y- m " < 0.4 GeV, (3.1) 
2m N ' v ; 

where rrijy denotes the nucleon mass and 

Py=Pn+Pj-P p (3.2) 

the 4-momentum of the undetected final state. In Eq. (3.2) p 7 and p p denote the 4-momenta 
of the incoming photon and the detected 7r + 7r~ pair and Pn is the 4-momentum of the struck 
nucleon which, for the calculation of py, is assumed to be at rest. In Fig. 2 we compare 
our calculation of the exclusive p° production cross section off hydrogen with experimental 
data [13-15]. In the whole Q 2 region covered by the HERMES experiment we find very 
good agreement for a broad range of the invariant mass W of the photon nucleon system. 
Also the slope of the differential production cross section is reproduced very well by the 
PYTHIA model as can be seen from the solid line in Fig. 3 where we show our calculation 
of ^{j*p -> p°p) together with the HERMES data [2]. 

For our calculations on exclusive p° electroproduction off nuclei we again use the kine- 
matic cuts of the HERMES collaboration [2,3]. This means that we restrict our exclusivity 
measure to the region 

-2 GeV < AE < 0.6 GeV, (3.3) 

and introduce a lower boundary for the four-momentum transfer \t — t max \ >0.09 GeV 2 as 
imposed by the HERMES collaboration to get rid of coherently produced p°. From the 
dashed line in Fig. 3 one sees that the differential p° electroproduction cross section off 14 N 
is again in excellent agreement with the HERMES data [2] . Throughout our calculations the 
effect of the nucleon potential turns out to be negligible. This is reasonable since the involved 
energies are much larger than the typical binding energies which are in the order of a few 
MeV. A combined effect of Fermi motion and Pauli blocking on the incoherent differential 
production cross section is visible at \t — t max \ < 0.1 GeV 2 . At \t — t max \ 0.05 GeV 2 the 
differential cross section is about 25% smaller than the calculation without Fermi motion 
and Pauli blocking (dotted curve). The reason for this reduction is that in the case that 
the bound nucleon moves towards the incoming photon the outgoing nucleon might be Pauli 
blocked for small momentum transfers. Note that this Pauli blocking just means that the 
whole nucleus absorbs the transferred momentum, i.e. this event contributes to the coherent 
production cross section which we do not consider here. 

In Fig. 4 we show the transparency ratio Ta for exclusive p° production as a function 
of the coherence length l p = q~ l . The solid line is the result that one gets if one uses Eq. 
(2.12) and accounts for two-body correlations by making the substitution [18] 
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n(b,Zi)(l -j (q c \zi ~ A)) 



(3.4) 



in (2.12) with q c = 0.78 GeV. This is necessary to avoid unphysical contributions from 
processes where Z{ ~ z which would contribute for small values of the coherence length l p . 
In Ref. [18] we showed that this Bessel function parameterization yields a good description 
of shadowing in photoabsorption. In Eq. (2.12) we use for the total p°N cross section 
<7p0 = 25 mb and for the elastic part a e \ = 3 mb. These two values correspond to the pN 
cross sections used within the transport model for the involved p° momenta. 

The result of the transport model is represented by the open squares. For each HERMES 
data point [3] we have made a separate calculation with the corresponding v and Q 2 . In the 
case of 14 N the Glauber and the transport calculation are in perfect agreement with each 
other and the experimental data. This demonstrates that, as we have discussed in Ref. [4], 
Glauber theory can be used for the FSI if the right kinematic constraints are applied. For 
comparison we also show the result of Huefner et al. [8] (dashed curve) as well as the result 
that one gets when using the approximate expression by Yennie [11] including two-body 
correlations (dotted curve). The somewhat larger transparency ratio of Huefner et al. arises 
from the different density distribution that the authors of Ref. [8] use which reduces the effect 
of FSI. The difference of the Yennie result (dotted curve) and our calculation arises mainly 
from the different cross section in the FSI (second replacement in (2.13)). The different 
phase factor leads to the change in the transparency ratio at small coherence lengths, since 



After applying all of the above cuts, nearly all of the detected p° stem from diffractive 
p° production for which the formation time is zero. The 14 N data seems to support the 
assumption that the time needed to put the preformed p° fluctuation on its mass shell and 
let the wave function evolve to that of a physical p° is small for the considered values of 
Q 2 . Furthermore, the photon energy is too low to yield a large enough 7 factor to make 
the formation length exceed the internucleon distance and make color transparency visible. 
This conclusion is at variance with that reached in Ref. [6]. 

We now turn to 84 Kr where we expect a stronger effect of the FSI. Unfortunately there is 
yet no data available to compare with. As can be seen from Fig. 4 the transport calculation 
for 84 Kr gives a slightly smaller transparency ratio than the Glauber calculation, especially 
at low values of the coherence length, i.e. small momenta of the produced p°. There are two 
reasons for this: About 10% of the difference arises from the fact that within the transport 
model the p° is allowed to decay into two pions. The probability that at least one of the 
pions interacts on its way out of the nucleus is about twice as large as that of the p°. The 
other reason is that in the Glauber calculation (2.12) only the inelastic part of the p°N cross 
section enters whereas the transport calculation contains the elastic part as well. Thus all 
elastic scattering events out of the experimentally imposed t-window are neglected in the 
Glauber description. It is because of this t-window that also elastic p°N scattering reduces 
the transport transparency ratio shown in Fig. 4. Both effects are more enhanced at lower 
energies and become negligible for the much smaller 14 N nucleus. 



Qv = I 



— 1 
v 
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IV. SUMMARY AND OUTLOOK 



We have developed a method to account for coherence length effects in incoherent elec- 
troproduction off nuclei which allows us to distinguish between the initial state interactions 
of the photon (shadowing) and the FSI of the reaction products. We have shown that our 
result is equivalent to the exact result of Glauber theory if one treats the FSI as a purely 
absorptive one. We have then performed a coupled channel treatment of the FSI within a 
semi-classical transport model and calculated the transparency ratio for exclusive incoherent 
p° photoproduction off 14 N and 84 Kr. The result for 14 N is in agreement with experimental 
data and with the Glauber prediction. The latter shows that in the case of 14 N Glauber 
theory is applicable after the kinematic cuts of the HERMES experiment are applied. Since 
we do not use a formation time for diffractively produced vector mesons we deduce that 
one cannot see an onset of color transparency in the nitrogen data. For the 84 Kr target 
no experimental data is available to compare with. However, we find deviations from the 
simple Glauber model because of the finite life time of the p° and elastic scattering out of 
the kinematically allowed |t|-region. These effects should be taken into account when eval- 
uating the 84 Kr data in search of color transparency. As discussed by Kopeliovich et al. [6] 
one might see an onset of color transparency when investigating the transparency ratio as 
a function of Q 2 for fixed coherence length [3]. Including a formation time for exclusive V 
production is not straightforward within our model since a change of the VN cross section 
during the vector meson formation will influence both the initial as well as the final state 
interactions. This is planned for future work, when hopefully also the 84 Kr data has become 
available. 
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was supported by DFG. 



In the following we show the equality of Eqs. (2.12) and (2.14). Since the real part of 
the VN scattering amplitude was neglected in the derivation of (2.14) we also set ay = 
in Equation (2.12) and end up with: 
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xjl-oy J dzin(b,Zi)cos[q v (zi- z)]exp -^°V J 



--a v I dz k n(b,z k ) 



dzin(b, Zi)e tqvZt exp 

oo 

a lN ^ VN [ d 2 b\ [ dzn(b, z )e-°v elT ^ 



--oy / dz k n(b,z k ) 



} 



ay \ dzin(b,zi) / dz 2 n(b, z 2 ) cos [qv(zi - z 2 )] 



x cxp 



Z2 



-a^ el / dz'n(b, z 1 ) - -a v \ dz'n{b, z 1 ) 



Z'Z 



a 



+- 



X 



V 



dzn(b, z) 



dzifi(b, Zi)e iqyZi exp 



z oo 

-\ov J dz'n(b,z') - J dz'n(b,z f ) 



}• (Al) 



In the second step we have renamed the integration variables and rewritten the integral 
limits. The first term of (Al) already equals that the first term in (2.14). The exponent in 
the second term of (Al) yields 

-a™ l T Z2 (b) - l -a v (T Z1 (b) - T Z2 (6)) = -\ « e < - <$) T Z2 (b) - \a v T Zl (b) (A2) 
and the exponent of the third term in (Al) is 

~\ (*vT Zi (b) - T z (b)) - \aT l T z {b) = - l -oyT Zi {b) + \^T z ib). (A3) 
We now further manipulate the last term of (Al): 

oo z 2 

— OO |— oo 

OO Z Z 

= °-^ J dz J d Zl J dz 2 n(b,z)n(b, Zl )n(b,z 2 )e a v T ^ cos [q v (z 1 - ^)] e -W r *i(S)+T. 2 ®). (A4) 

— OO — OO — oo 

Rewriting the integral in the following form 
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OO 2 



dz / dz\ l dz. 



-OO — OO — OO 



OO 22 OO OO 21 OO 

cb 2 J dz\ J dz + y d^i y cb 2 y dz 

OO — OO 22 — OO — OO 21 



(A5) 



and using the symmetry of the integrand with respect to the variables z x and z 2 we can 
rewrite (A4) as 



^ y dz2n(b,z 2 ) y ^in(Mi)cos[g K (^i-^2)]e-^ T «i® +T *2( s )) y dzn(f,z)e^ 

— OO — OO 22 

OO 2 2 

= ~\^ J dz Mb, *i) j d Zl n(b, z 2 ) cos [q v ( Zl - z 2 )\ e -^ T ^+ T ^)) (i _ e ^ 2 ^ 

— OO — OO 

OO OO 

= -)~ l I dz 1 n(b,z 1 ) Jdz 2 n(b,z 2 )cos[q v (z 1 -z 2 )]e-^ T ^ +T ^ 

— OO 21 
OO OO 

+7^ / <km(Mi) f dz 2 n{b,z 2 )cos[q v { Zl - ^ 2 )] e^" 6 ~ CT vW&V V^iW). (A6) 

2 °V J J 



Zl 



In the first step we have performed the integral over z. In the last step we have again 
rewritten the integral over z\ and z 2 to be able to combine the second term of (A6) with 
the second term of (Al) to get the second term in (2.14). Now we only have to show the 
equality of the last term in (2.14) and the first one in (A6). Therefore we rewrite 



A (4 



1_<± 



dzn(b,z)e iqvZ e- 1 2 avT *® 

OO 

dz x [ dz 2 n(b,z 1 )n(b,z 2 )cos[q v (z 1 - ^)] e -W( T «i®+ T « 2 ®) 



(AT) 



-oo — OO 



and reformulate the integral as 



OO OO 

y dz\ y dz 2 



OO OO OO OO 



dz\ \ dz 2 + 



I dZ2 I 



dz\ 



-OO — OO — OO 21 — OO 22 

OO OO 



2 J dz\ J dz 2 

— OO 21 



(A8) 



where the second equality again follows from the symmetry of the integrand with respect to 
the exchange of Z\ and z 2 . The two expressions (2.12) and (2.14) are therefore equivalent. 
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FIGURES 
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FIG. 1. PYTHIA result for the fraction of the total j*N cross section that is contributed by 
the VMD (solid line), GVMD (dashed) and direct (dotted) part of the photon as a function of Q 2 
for fixed invariant mass W = 5 GeV. Note that in contrast to the authors of Ref. [9] we call every 
pointlike interaction of the photon a direct process. 
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FIG. 2. The virtual-photoproduction cross section for p° production plotted versus the invariant 
mass W at average Q 2 values of 0.83, 1.3, 2.3 and 4.0 GeV 2 . The data has been taken from the 
HERMES [13] (filled circles), E665 [14] (open circles) and CHIO collaboration [15] (open diamonds). 
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FIG. 3. The exclusive incoherent differential p° production cross section for l H (solid line) and 
1 N (dashed line) in comparison with experimental data from the HERMES collaboration [2]. The 
calculation has been performed for v = 13 GeV and Q 2 = 1.7 GeV 2 . The dotted curve represents 
a calculation for 14 N without Pauli blocking and Fermi motion. For details on the exclusivity 
measure see text. 
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FIG. 4. Nuclear transparency ratio T4 for p° electroproduction plotted versus the coherence 
length of the p° component of the photon. The data is taken from [3] (filled circles), [16] (open 
circles) and [17] (open diamonds). The solid line represents the Glauber result when using (2.12) 
including nucleon-nucleon correlations, the dotted line shows the effect of the substitutions (2.13). 
The dashed curve represents the Glauber result using (2.14) with the density distribution of Ref. [8]. 
For each transparency ratio calculated within our transport model (open squares) we used the 
average value of Q 2 and v of the corresponding HERMES data point. 
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